Magnetic order and ice rules in the multiferroic spinel FeV204. 
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We present a neutron diffraction study of FeV204,, which is rare in exhibiting spin and orbital 
degrees of freedom on both cation sublattices of the spinel structure. Our data confirm the existence 
of three structural phase transitions previously identified with x-ray powder diffraction, and reveal 
that the lower two transitions are associated with sequential coUinear and canted ferrimagnetic 
transitions involving both cation sites. Through consideration of local crystal and spin symmetry, 
we further conclude that Fe'^'^ cations are ferro-orbitally ordered below 135K and V^^ orbitals 
order at 60K in accordance with predictions for vanadium spinels with large trigonal distortions and 
strong spin-orbit coupling. Intriguingly, the direction of ordered vanadium spins at low temperature 
obey 'ice rules' more commonly associated with the frustrated rare-earth pyrochlore systems. 
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Magnetic insulators with orbital degrces-of-freedom 
are of great interest to both the condensed matter and 
applied physics communities. Transition-metal spinels 
are a particular example where orbital degeneracy and 
the concomitant orbital ordering have been shown to 
profoundly affect both th e cr ystal structure and mag- 
netic exchange interactional^. In recent years, the ma- 
terial FeV204 has emerged as an important model sys- 
tem to explore this physics. A rare example of a spinel 
with two orbitally- active cation sites, FeV204 is a mul- 
tiferroic, cont aining both ferroelectric!^ and ferrimag- 
netic momentsDSIMl at low temperatures. Moreover, 
separate experiments have demonstrated that the ap- 
plied magnetic fields can both suppress the ferroelectric 
moments^ and induce structural detwinnin^. Previous 
x-ray diffraction measurement^ mapped out a sequence 
of structural phase transitions, however despite the obvi- 
ous importance of magneto-elastic and magneto-electric 
coupling, magnetic structures have remained unknown. 
In this Letter, we present new neutron scattering data 
definitively solving the magnetic structures as a func- 
tion of temperature. The results establish that FeV2 04 
contains unexpected and very interesting parallels to the 
well-known spin-ice problem of frustrated magnetism^, 
realized at low temperatures in selected rare-earth con- 
taining pyrochlore materials. 

The spinel vanadates, AV2O4, with divalent A-site 
cations have been central to the study of orbital order 
and the role of orbital degrees-of-freedom in geometri- 
cally frustrated antiferromagnets. In these systems, the 
pyrochlore sublattice of the spinel structure is occupied 
by 1/^+ (3d^) cations in an octahedral oxygen environ- 
ment, giving spin S=l and an orbital triplet degree-of- 
freedom. Materials such as ZnV20jHini, CdV204" and 
MgY20p^, where the A-site is magnetically neutral, 
are typified by two transitions: an orbital ordering transi- 
tion which lowers crystal symmetry at an upper tempera- 
ture, followed by a second magnetic transition to a state 



with a complex pattern of ordered spins. In MnV204, 
with spin-only Mn^+ on the A-site, the system first or- 
ders magnetically to a collinear spin state involving both 
cation sites, only to develop a non-coUinear structure 
below a second transition at a lower temperatur^^JtiS] 
Despite being electronically and structurally similar at 
high temperatures, the low temperature properties of 
these materials are strikingly different, as the relative im- 
portance of spin-orbit coupling, tetragonal and trigonal 
crystal fields vary in different systems. In each of these 
materials, debate continues and centers on understand- 
ing the nature of the orbital order inferred from mea- 
sured crystal symmetry and magnetic structur^i^HlIl^ 
largely separate line of research has concentrated on the 
role of orbitals on the A-site sublattice, for example in 
FeCr204'2SH2Slor FeSc2Sj^213iI however, very few materi- 
als have been identified with orbitally active ions on both 
cation sites of the spinel structure. 

In this regard, FeV2 04 is an extraordinary material. In 
addition to V^^ on the pyrochlore sublattice, this system 
contains tetrahedrally-coordinated Fe^^ (3d^) cations on 
the A-site of the spinel structure with high-spin S=2 and 
two-fold orbital degeneracy. As might be expected, fur- 
ther interactions result in additional complex and inter- 
esting behavior. Specific heat measurements have re- 
vealed three phase transitions in this material with 
= 138, 111 and 56 RH. Powder x-ray scattering has 
associated these temperatures with a cascade of struc- 
tural transitions, from cubic to high-temperature tetrag- 
onal (HTT) to face-centered orthorhombic (FCO) to low- 
temperature tetragonal (LTT) with decreasing tempera- 
ture. Both Mossbauei'^ and magnetizatiorP^l measure- 
ments indicate that at least one of these transitions 
(HTT— >FCO) additionally involves the ordering of mag- 
netic degrees-of-freedom. A fourth structural phase tran- 
sition to low-temperature orthorhombic (LTO) symmetry 
has been reported in one single crystalline studjpl, and at 
least one AC-susceptibility study has identified a glass- 
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like transition at 85 kP. Despite these advances and the 
implied coupling of electric, lattice and magnetic degrees- 
of-freedom, very little is known about the local spin and 
orbital structure in the distinct phases of FeV204. 

In this Letter, we present results of a neutron pow- 
der diffraction (NPD) study of FeV204 at tempera- 
tures below 200 K. Wc confirm the cubic-HTT-FCO- 
LTT sequence of structural phase transitions reported 
by powder x-ray scattering, but see no evidence for 
the purported fourth transition to a low-temperature or- 
thorhombic phase- . Hitherto unreported, we discuss the 
evolution of the crystal field environment of V^'^ cations, 
and in particular present evidence for a large trigonal dis- 
tortion of the octahedral oxygen cage about these sites. 
We additionally report the existence of two separate mag- 
netic ordering transitions: a collinear spin-ordering co- 
incident with the HTT-FCO transition and a tilting of 
V^^ moments to point along the <111> crystal direc- 
tions when symmetry changes from FCO to LTT. We dis- 
cuss the implications of the current observations for the 
orbital order, and in particular argue that this material 
is well-described by models which assume large trigonal 
distortion and strong spin-orbit coupling. 

Poly crystalline samples of FeV204 were prepared by 
solid state reaction, and NPD measurements were per- 
formed using the HB2A diffractometer at the High Flux 
Isotope Reactor!^. Further elastic neutron scattering data 
was taken on a crystal grown via the float zone method 
using the CG-4C Cold Neutron Triple- Axis Spectrometer 
of HFIR. Details about sample preparation and experi- 
mental method can be found in the Supplementary Ma- 
terial, along with characterization data on both powder 
and single-crystal samples. 

The results of our high-resolution neutron diffraction 
measurements at 200 K, 125 K, 75 K and 4 K are shown 
in Fig.[3]for 20 < 85°. The highest temperature data con- 
firm the spinel structure at these temperatures and the 
near phase purity. Refinement implies that a fraction 
of the vanadium sites is occupied by iron cations, giving 
chemical formula FejFeo.ossVo. 942)204 ■ This is similar to 
other samples in the literature^. There is also a small 
temperature-independent impurity phase, possibly asso- 
ciated with an oxide of vanadium. The temperatures 
of remaining high-resolution patterns were chosen to be 
well within each of the five phases identified by Katsu- 
fuji et al. with single-crystal x-ray diffractiorFi. Con- 
trary to those results however, no discernable difference 
was seen between the patterns at 50 K and 4K, and so 
only the 4K pattern is plotted here. The four unique 
patterns were well-described by the space groups Fdim, 
Ml/amd, Fddd and 141 /amd with decreasing temper- 
ature, consistent with previous analysis and reflecting 
the cubic-HTT-FCO-LTT sequence of structural phase 
transitions^. Patterns to 29 = 154° and full refinement 
information are given as Supplementary Material, along 
with a more in-depth discussion of the absent LTT-LTO 
phase transition. 

The sequence of observed structural transitions is il- 
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FIG. 1: Plots of raw NPD data (circles) measured at T = 
200 K, 125 K, 75 K and 4 K for 261 < 85°. Solid lines are 
results of Rietveld refinements described in the main text. 
Differences between observed and calculated intensities are 
shown directly below the respective patterns. See the main 
text and Figure |2] for more detail. 



lustrated more clearly in Fig. [ga) and (b). Fig. [ga) 
shows the successive splittings and recombination of the 
cubic (004) peak for each of the four patterns. Fig. [2|b) 
shows a plot of the lattice parameters extracted from 
lower-resolution NPD data assuming the crystal symme- 
tries determined above. (For the sake of comparison, all 
lattice parameters plotted here are those for the cubic 
unit cell.) Transition temperatures are identified as 135 
K, 110 K and 60 K, within an uncertainty of 5 K. These 
estimates are similar to the values inferred from heat ca- 
pacity in recent measurements of powders by Zhang et 
al^^, as well as our own single crystal data (see Supple- 
mentary Material). 

Fig. ^c) and (d) demonstrate the neutron scatter- 
ing intensity of the cubic (111) peak as a function of 
temperature. Contrary to what is expected from struc- 
tural considerations alone, NPD reveals a significant in- 
crease in this peak below the HTT^-FCO transition and 
again below the FCO^LTT transition. At both tem- 
peratures, an increase in net moment has been noted 
in bulk magnetization datcP^^I. For comparison, we also 
show in Fig. ^d) the temperature evolution of the (220) 
peak from NPD and the (much weaker) (200) peak from 
single-crystal elastic scattering. The (220) peak carries 
only information about the iron sublattice, and its tem- 
perature dependence refiects only the 110 K transition. 
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FIG. 2: (a) An expanded view of the NPD data around the 
cubic (004) Bragg position for each of the patterns in Fig. [s] 
(b) Plot of the cubic lattice parameters extracted from NPD 
data, (c) An expanded view around the cubic (111) Bragg 
position for the same temperatures as (a), (d) A comparison 
of the temperature evolution of NPD (111) and (220) Bragg 
peak intensities, above values their values at 180 K. Overlaid 
is the scaled intensity of the (200) Bragg peak, as measured 
in a single crystal sample. 



Conversely, the intensity of the (200) peak is zero at high- 
temperature and grows below the transition at 60 K. 
These observations are consistent with our Rietveld re- 
finements, which associate two lowest phase transitions 
with the onset of collinear and canted ferrimagnetism, as 
discussed in more detail below. 

To gain insight about the crystal field environment at 
the Fe^+ and V'^^ sites, we further plot in Fig. [s] the 
O-cation-O angles as a function of temperature for the 
FeOi tetrahedra and VOq octahedra, respectively. In 
the high-temperature cubic (Fd3m) phase, the interior 
angles of the FeO^ tetrahedra (Fig. [Sjja)) are sixfold- 
degenerate with a value of 109.4°. This reflects the fact 
that the tetrahedra are regular, and thus the orbital 
ground-state of Fe'^'^ is two-fold degenerate. At the 135 
K transition, this degeneracy is lifted by a uniform com- 
pression in the cubic (001) direction, lessening four of 
the six O — Fe — O angles and lowering the energy of the 
single-electron (i^2 orbital relative to the d^'^-y'^ orbital. 
As noted previously by Katsufuji et al^ and Sarkaret 
al^, this behavior can be associated with ferro-orbital 
order on the Fe^'^ site. The orthorhonibic transition in- 
troduces an additional in-plane compression of the tetra- 
hedra, which increases with decreasing temperature un- 
til the system locks-in to the low-temperature tetragonal 
structure. Notably, these subsequent transitions affect all 
Fe^"*" equivalently and are consistent with ferro-orbital 
order persisting to lowest temperatures (though the exact 
distribution of electron density will evolve upon cooling) . 
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FIG. 3: Plots versus temperature of the interior O — M — O 
angles (M £ {Fe^+ ,V^+)) of the (a) FcOa tetrahedra and 
(b) VOq octahedra, shown schematically at the top. In pic- 
tures, bonds of equal length are shown with equivalent color- 
ing. Arrows denote the largest distortion of each polyhedron: 
a compression of the FeOi tetrahedron at 135 K along [001] 
and a trigonal distortion of the VO^ octahedron along <111> 
directions that exists at all temperatures. 



Conversely, in the Fdim phase the octahedra dis- 
play a significant trigonal distortion from cubic symme- 
try. This is evident in Fig. |3][b), where one can see that 
at high temperatures the twelve interior angles of the 
VOq octahedra are split into two inequivalent branches of 
six with values significantly different from 90° (expected 
for regular octahedra). This splitting corresponds to a 
stretching of the octahedra in alternating <111> direc- 
tions along chains in the [110] direction. The distortion 
can be associated directly with the oxygen position pa- 
rameter within the FdZm space group. The angular sep- 
aration at high temperatures (10.7°) is comparable to 
what is inferred from the structures of MnV204 (12.9°) 
and MgV204(8.7°), where the importa nce of tr igonal dis- 
tortion has already been acknowledgecP^'^ilSSi^ Katsufuji 
et al. previously noted a tetragonal compression of the 
VOq octahedra at the 135 K transitiorP. This is observ- 
able in Fig. [3]jb) as well. Notably however, the trigonal 
splitting of the 0-V-O angles in the cubic state is nearly 
ten times larger than the perturbations at lower temper- 
atures. (Note the break in the vertical axis.) Thus our 
data imply that the trigonal distortion is the dominant 
non-cubic term in the crystal field environment of V'^'^ . 

Information about the determined magnetic structure 
is displayed in Fig.|4] As h ypo thesized previously in anal- 
ogy to other spinel systemJI'S]^ the HTT— J-FCO structural 
phase transition coincides with the onset of collinear fer- 
rimagnetism, wherein Fe^^ moments align ferromagnet- 
ically along the longest orthorhombic axis and mo- 
ments align antiparallel. This is the ordered state fa- 
vored by antiferromagnetic superexchange between the 
Fe^^ and V^^ cations and implies that the transition is 
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driven by such magnetic interactions. The orthorhombic 
structural transition then hkely results from a coupling 
between cation spins and phonon modes which strongly 
affect superexchange paths, similar to what has been sug- 
gested by recent optical measurements of FeCr2 04 and 
NiCraOjSal. As shown in Fig. ^ a), the ordered moment 
sizes on both sublattices grow below T/vi = 110 K, with 
the vanadium moment perhaps seeing additional increase 
at 60 K. The ordered moment sizes at lowest tempera- 
tures are 4.0/xb and 0.85^_b on the Fe^"*" and sites, 
respectively. Whereas the number for Fe^'^ is in line 
with expectations for a S=2 system, the value for V^'^ 
is significantly less than the 2.0 hb expected for an S=l 
cation with quenched orbital angular momentum. 

Refinements indicate no tendency for V^'^ to cant away 
from the axis defined by the i^e^"*" moments in the or- 
thorhombic phase, though symmetry would allow such 
behavior. Nonetheless, it is shown in Fig. |4][b) that the 
refined angle between the V^^ moment and the Fe'^^ 
spin axis jumps abruptly from zero to ^ 55°±4° below 
the FCO-LTT transition at Tn2 = 60 K. This canting 
transition is illustrated in Figs. |4]Jc) and (d), where we 
sketch the ordered states above and below Tn2 as viewed 
along the cubic (100) direction. The existence of a cant- 
ing transition is further supported by the single-crystal 
data in Fig. ^d) , as the intensity of the (200) Bragg peak 
is proportional to the square of the in-plane vanadium 
moment in this model. The low-temperature value of 
55° for the canting differs from the 33° inferred recently 
by Kang et al. from NMR data'^^'. However as noted 
previousljEHj a value of 33° cannot be reconciled with the 
large ferrimagnetic moment observed via magnetizatiorP 
(which implies an angle closer to 60°). We feel that the 
neutron diffraction data provide the more reliable esti- 
mate. 

The evolution from a collinear to canted ferrimagnetic 
state is reminiscent of the magnetic behavior of MnV204, 
where the canting was accompanied by the gapping of 
magnetic excitationJi^. Significantly though, we note 
that the canting angle observed here is smaller than the 
65° reported for MnV20j^ and is equal within error 
to the 54.7° between the (001) and (111) directions in 
a cubic lattice. In addition to the canting angle, the 
'in-plane' (i.e. _L(001)) spin directions imposed by the 
lAl/amd space group are also different from what has 
been assumed for the ordered state of MnV20ii3, and 
in fact refinement suggests the low-temperature phase 
of FeV204 is characterized by a '2-in-2-out' structure 
on the pyrochlore lattice (see Fig. |4] (e)), more com- 
monly associated with 'ice-rules' in frustrated rare-earth 
antiferromagnetJ^. 

In MnV204, the low temperature canting transition 
has been associat ed with the ordering of V^^ orbital 
degrees-of-freedorrfi^l^^tlll, and similar physics is likely 
at play here. Theoretically, it has been shown that the 
ground state properties of both MnV20i23and FeV20j2S 
are strongly infiuenced by the alternating pattern of trig- 
onal distortion axes which now has been observed in 
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FIG. 4: (a) Plot of the ordered moment size on the Fe^+ 
and V'^'^ cation sites and (b) the angle between the V^^ mo- 
ments and the axis defined by Fe^^ moments, (c) A sketch 
of the collinear ferrimagnetic state in the FCO phase, as seen 
along the cubic (001) direction, and (d) a similar sketch of the 
canted state seen in the LTT phase. Details about the order- 
ing patterns are given in the main text, (e) A view of V'^^ mo- 
ments in low-temperature ordered state clearly demonstrating 
the 2-in-2-out spin structure on the pyrochlore sublattice. 



both materials. The resultant physics is captured by the 
'quantum 120° model' analysis of Chern et aP^, who con- 
sidered vanadium spinels with dominant trigonal distor- 
tions and showed for these materials that an orbitally or- 
dered state is stabilized at low temperature with electron 
density modulated along <110> directions. The spin or- 
der is then determined by a competition between rela- 
tivistic spin-orbit coupling and orbital exchange between 
vanadium sites. When the latter wins out, the ordered 
spins tend to lie along primary cubic axes in the plane 
perpendicular to the ordered A-site spins. When spin- 
orbit coupling is dominant, a strong single- ion anisotropy 
forces the spins to lie along cubic <111> directions and 
favors the '2-in-2-out' structure seen in Fig.[4j It has been 
shown that MnV204 is well-described by an intermediate 
case, where spin-orbit coupling and orbital exchange are 
comparabl^^H, xhe present data suggests that FeV204 
is well within the strong spin-orbit coupling limit of this 
model. It is worth noting that strong spin-orbit coupling 
would have the additional effect of reducing the V^^ or- 
dered moment size from spin-only values, consistent with 
observation. 

From a pure physics perspective, the '2-in-2-out' spin 
structure itself may be one of the most interesting new 
observations. In the rare-earth pyrochlore antiferromag- 
nets, a local '2-in-2-out' constraint arises from a net ferro- 
magnetic coupling between Ising spins with quantization 
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axis along <111> and leads to the macroscopic ground- 
state degeneracy underlying the much-studied 'spin-ice' 
class of compoundJ^. Our data imply that in FeV2 04 
an equivalent Ising condition on vanadium spins results 
from the large trigonal distortion of the VOq octahedra 
in the presence of strong spin-orbit coupling. Thus, one 
way of interpreting the low-temperature magnetic state 
in this material is as a c?-electron analogue of the spin- 
ice problem, but persisting to temperatures as high as 
60 K. In such a picture, the role of the large ordered 
spin on the Fe'^^ site is to act as an applied field and 
break the macroscopic degeneracy; indeed, the magnetic 
order in the LTT phase is exactly what has been ob- 
served in spin-ice materials Ho^Ti^O'j and Dy2Ti20-! 
with field applied along (001^^. An analogous situa- 
tion may be found in the 'ordered spin-ice' phases of 



SmaMoa 0,1221 and NdaMoa 0,321, where rare-earth mo- 
ments on the pyrochlore lattice couple to ferromagnet- 
ically ordered Mo spins. In Tb2Sn207, a similar 'or- 
dered spin-ice' state has also been seerP^]^ and was ini- 
tially associated with ferromagnetic near-neighbor cou- 
pling in the presence of moderate single-ion anisotropy. 
The results reported here provide an important link be- 
tween such rare-earth pyrochlore antiferromagnets and 
trigonally-distorted spinels such as FeV204. 
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SUPPORTING MATERIALS 



I. SAMPLE PREPARATION AND 
CHARACTERIZATION 

The samples investigated in this manuscript were pre- 
pared by grinding appropriate mixtures of Fe2 03, Fe and 
V2O3, pressing under 400 atm hydrostatic pressure and 
calcining in flowing argon at 950 C for two days, with in- 
termediate grinding to ensure sample homogeneity. The 
sample was characterized by bulk magnetization using a 
Quantum Design Magnetic Property Measurement Sys- 
tem (MPMS), and the results are shown in Figure SI. Of 
note are the distinct anomalies at the T = 110 K and 60 
K, consistent with the ordered states inferred from our 
neutron scattering data. 

A single crystal of FeV204 was also recently been pre- 
pared via the float-zone method. The crystal growth 
was carried out in argon in an IR-heated image furnace 
equipped with two halogen lamps and double ellipsoidal 
mirrors. Feed and seed rods rotated in opposite direc- 
tions at 25 rpm during crystal growth at a rate of 25 
mm/hr. In addition to magnetic susceptibility, a por- 
tion of the crystal was characterized via heat capacity 
using a Quantum Design Physical Property Measurement 
System (PPMS). The results of these measurements are 
shown in (Fig. S2) indicate second-order transitions with 
Tc = 138, 107 and 60 K, consistent with our neutron pow- 
der diffraction results and heat capacity measurements on 
powder^. Notably, there is no sign of a transition at 35 
K, as reported by single crystal x-ray scattering studiefP. 
This is discussed in more detail below. 



II. NEUTRON SCATTERING 
MEASUREMENTS 

Powders were investigated using neutron powder 
diffraction, using the HB2A diffractometer at the High 
Flux Isotope Reactoi^. High resolution data were col- 
lected at 200 K, 125 K, 75 K, 50 K and 4K using a 12'- 
21'-6' collimation and two different wavelengths (1.538 
A and 2.41 A). Additional lower resolution measure- 
ments were carried out using the wavelength 1.538 A and 
12'- 40'-6' collimation to monitor the change in the scat- 
tering across various phase transitions. Rietveld refine- 
ments were performed on combined data sets from two 
separate wavelengths using FULLPROlS. Data over a 
limited range of 29 and analysis are presented in the main 
text of this Letter. Diffraction patterns covering the en- 
tire accessible range of 29 are shown in Figure S3, and 
structural parameters inferred from Rietveld refinements 
are given in Table [l] 

The single crystal was investigated using elastic neu- 
tron scattering. Measurements were performed using the 
cold neutron triple-axis spectrometer (CTAX) of HFIR 
with guide-open-80'-open collimation. The energy of the 



scattered neutrons was fixed at = 5 meV. Higher order 
contamination was removed by a cooled Be filter placed 
between the sample and analyzer. Elastic scattering data 
pertinent to the current discussion were presented in Fig- 
ure 2 of the main text. Further data will be presented in 
a future publication. 



III. ABSENCE OF A LOW-TEMPERATURE 
ORTHORHOMBIC PHASE 



Figure S3 (and Figure 1 in the main text) shows high- 
resolution neutron scattering patterns for temperatures 
200, 125, 75 and 4 K, chosen to be within each of the 
distinct phases identified by the x-ray scattering study 
of Katsufuji et al? . As mentioned in the main text, a 
high-resolution pattern also exists for T = 50 K, but 
was indistinguishable from the pattern at T = 4 K. This 
is consistent with Figure S2 above and Figure 2(a) of 
the main text, which show sign of a transition below the 
FTO— J-LTT transition at 60K. This is in apparent con- 
tradiction, however, to the paper of Katsufuji et al., who 
identified a fourth, LTT— >-LTO structural phase transi- 
tion in a single crystal with Tc = 35 M^. With regard to 
this point, it is interesting to note no sign of a LTT— J^LTO 
transition was observed in either the powder diffraction 
data of Katsufuji et al^ or the heat capacity data of 
Zhang et al^. We further note that absence of a peak at 
35 K in our single crystal heat capacity data shown in 
Figure S2. This last observation and the magnitude of 
the peak splitting reported in Ref [2] seem to preclude the 
idea that the LTT— t'LTO transition is merely a subtle 
effect which is detectable only in single crystal samples. 

An alternative explanation is that small differences in 
the composition of different samples could be responsible 
for the discrepancy in low-temperature behavior. It has 
been noted on several occasions that samples of FeV204 
display a stro ng te ndency towards excess iron on the 
vanadium siteJ^Ml (a,n effect also seen by us), and de- 
tailed study by Rogers et al. has demonstrated that even 
small variations in iron stoichiometry can significantly af- 
fect the physical properties of this systerrPi. This includes 
a change in the number of phase transitions observable 
with x-ray scattering above nitrogen temperatures (the 
limit of that particular study). Though a similar study 
has not been done at low temperatures to our knowl- 
edge, this material-specific sensitivity to iron concentra- 
tion makes sample stoichiometry a prime candidate to 
explain the variation in the presence of the LTT->LTO 
transition, even within Ref. [2] We point to the variation 
in FCO-LTT transition temperatures between the pow- 
der (Tc = 60 K) and crystal (T^ = 85 K) samples of Ref.H 
as further support for this idea. 
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IV. GLASS-LIKE BEHAVIOR IN THE FTO 
PHASE 

We also noted briefly in the main text that a spin-glass 
transition has been identified by Nishihara et al. at 85 
K using AC-susceptibilitji^. Glassy correlations would 
appear in neutron measurements as diffuse scattering in- 
tensity centered around the locations of specific Bragg 
peaks. In our powder measurements of FeV204, we did 
observe such scattering about the (111) position, however 
closer inspection reveals that the diffuse scattering inten- 
sity peaks at the 110 K tetragonal-orthorhombic phase 
transition, decreasing at lower temperatures with no fur- 
ther signature at 85 K. This is shown below in Figure S4. 



The more obvious interpetation of our diffuse scatter- 
ing then is as critical correlations associated with a fer- 
rimagnetic phase transition in a frustrated system, and 
not as a signature of glassiness. We do note however 
that these observations do not preclude the possibility 
of weak glassy dynamics associated with a different re- 
ciprocal lattice position, such as (200), but the intensity 
of the (111) scattering here makes powder diffraction an 
unsuitable vehicle for exploring such physics. In fact, the 
single crystal data presented in Figure 2(d) of the main 
text reveals significant correlated paramagnetic behavior 
above the 60 K transition which might be associated with 
the 'glassy' dynamics of Rcf. El 
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FIG. 1: Zero-field cooled (ZFC) and field-cooled (FC) measurements of bulk magnetization on powders of FeV204. Details are 
in text of Supplementary Information. 
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FIG. 2: Heat capacity of single crystalline FeV204, as described in Supplementary Information text. Peaks imply transitions 
at 138 K, 107 K and 60 K, consistent with our results in the main text. 
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FIG. 3: Plots of raw NPD data (circles) measured at T = 200 K, 125 K, 75 K and 4 K. Solid lines are results of Rietveld 
refinements described in the main text. Differences between observed and calculated intensities are shown directly below the 
respective patterns. 
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TABLE I: Structural parameters extracted from simultaneous Rietveld refinement of A = 1.538A and A = 2A1A NPD data. 
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FIG. 4: (left) A plot of the weak diffuse scattering intensity about the (111) Bragg position for each of the patterns shown 
in Figure S3, (right) The fit intensity of the diffuse scattering, as compared to the (111) Bragg intensity around the 110 K 
ferrimagnetic transition. 



